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Abstract 

We present a statistical model of a meandering river on an alluvial plane 
which is motivated by the physical non-linear dynamics of the river channel 
migration and by describing heterogeneity of the terrain by noise. We study 
the dynamics analytically and numerically. The motion of the river channel 
is unstable and we show that by inclusion of the formation of ox-bow lakes, 
the system may be stabilised. We then calculate the steady state and show 
that it is in agreement with simulations and measurements of field data. 
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Meandering rivers are ubiquitous in nature and can be found all over the world. It is a 
signature of their intellectual challenge that they have captured the imagination of scientists 
in quite differing fields [|I|. They are important not only as examples of a system far from 
equilibrium undergoing an interesting dynamics but because abandoned channels of the main 
river become silted up and over thousands of years these silted channels become compressed 
to form shale beds between which oil is often found. Previous models that have tried to 
study river meanders fall into two broad classes. The first could be considered as empiric 
and are justified by the fact that Leopold and Wolman showed by an analysis of field 
data from around the world, that there was a nearly constant ratio of radius of curvature to 
meander length and of radius of curvature to meander width. This 'explains' the visual 'scale 
invariance ' when one observes rivers. They therefore concluded that " meander geometry is 
related in some unknown manner to a more general mechanical theory" . Consequently there 
have since been a number of papers that modelled meanders as variants of random 
walks with essentially no justification. The second type is characterised by an analysis of 
the complex flow inside the channel and as a result yields relatively little information 
about large scale or long time features of the river pattern. In this letter, a slightly different 
approach is taken; we attempt to find the simplest possible model of the behaviour of a 
meandering river which still retains the essential physics. The model is statistical in the 
sense that there is a spatio-temporal random aspect to the dynamics. Statistical models 



have been introduced in the study of river networks [p!0|-[T2[| but have not been used in 
the study of single channels. Also, to the authors knowledge, a stable state has not been 
constructed for any dynamic models of river meanders. 

We consider the river on the length scale such that its channel may be considered to be 
a curve of length L where L —>■ oo, in two dimensions, r{s, t) = [x{s, t), y{s, t)] where s refers 
to the arc length position on the curve. We include all the major features of the dynamics 
of the channel and using them calculate the steady state. The derivation is outlined here 
and details will be reported elsewhere []T^. Our theory of meanders therefore acts as a 
connection between both previous approaches. 



To calculate the equation of motion (EOM) of the river 'curve', we must begin by an 
analysis of the flow of water in a curved channel. The channel is considered to be narrow 
with width b and we can describe the coordinates of the channel as s and n corresponding 
to moving along the centre-line of the channel and perpendicular to it, and note n = n{s). 
We transform the Navier-Stokes equation for incompressible flow to the curved channel 
coordinates system. Using the similarly transformed continuity equation and assuming no 
stress at the top surface and bed stresses at the bottom of the channel, it can be shown 
using simple arguments |^ that the major effect of this is that a velocity gradient develops 
across the channel with the velocity of water at the outer bank being greatest. That is if the 
outer bank corresponds to n = 6/2 and the inner bank n = —b/2, then the velocity of fluid 
in the channel is given to leading order by u{n) ^ ^0(1 + ^ i^i^)) where k,{s) = \d'^r/ds'^\ is 
the local curvature and Uq is a constant. We do not include gravity so the plane is flat but 
we do have a pressure gradient driving the flow through the channel. In the analysis above 
the approximation is made that the channel is locally smooth, i.e there are no sharp changes 
in direction meaning the local curvature is small. Assuming erosion on the outer bank 
is balanced by deposition on the inner bank and that the rate of erosion on the outer bank 
is proportional to the velocity, the rate of normal, i.e. perpendicular to the curve, migration 
of the channel can be described by the relationship 

Rl cx Kis) . (1) 

The condition that the curve is locally smooth is satisfled by adding another piece to the 
rate of migration; i.e. 



+ (2) 



This is obtained from the Rayleighian dynamics [|I^ of a 'global' free energy which minimises 
curvature along the whole meandering river curve, jF[r] = / ds\d'^r/ds'^\'^. Evolution of 
curves undergoing similar dynamics as these flrst two bits of Rm have been studied in the 
context of crystal growth ||15|,|16| and chemical front motion [jl^. Finally, we have a dynamic 



noise term to model heterogeneities in space, such as different types of rock, and in time, 
for example, variation in flow of the river and re-eroding parts of the plain that have been 
eroded previously. 



Dynamic noise is more appropriate than quenched noise because the river is constantly 
reworking its plain as it meanders and over long periods changing the nature of the terrain. 
For simplicity we take Gaussian uncorrelated noise. 

To get the correct EOM we must appeal to reparametrisation invariance because the 
arc length, s is not constant in time. Physically, it is known that the normal velocity of a 
curve, in this case n ■ dtr{s,t) = I]i=o'^i-^m 5 is the unit normal) is independent of its 
parametrisation so for simplicity the EOM may be considered flrst in the 'normal gauge' 
where the tangential velocity of the curve is zero. Then using the Frenet formulae from 
differential geometry the EOM in the 'arc length gauge' may then be calculated to give. 



where the normal velocity may be written V = Ak,{s) — B (k^ — d'^n/ds^) + ri{s,t) = — n • 



{Ad'^r/ds'^ + Bd^T / ds^)+rj[s, t) . The smallest length-scale of the system is therefore \J B / A. 

Recalling that dr/ds = s, the unit tangent vector to the curve, one sees that the second term 
on the Ihs of equation(^) is a non-local tangential velocity corresponding to a 'stretching' 
or lengthening of the curve due the normal growth. To calculate the effective motion of 
the curve we replace the tangential velocity by its average, i.e. because we are dealing 
with a Langevin equation with the prescence of a noise term, and we assume the system 
will always be near some steady state, we make the pre-averaging approximation where 
complicated parameters in the equation are replaced by their mean values , so we define 
C{A, B) = nVds') = J Vv KVds'P{[r],t), and consequently we may write the EOM of 



(3) 
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the river as 
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where 

{r]{s, t))=0, {r]{s, t)r]{s', t')) = D6{s - s')6{t - t') (6) 

We have thus calculated the first part of the problem, the local dynamics of the meander- 
ing river as a non-linear noisy equation. Noisy non-linear equations have been widely studied 
in the context of surface deposition |[TB| particularly when they lead to critical behaviour. 
Unlike these models meandering river systems do evolve into statistically invariant stable 
states. This equation also bears some superficial similarity to the Rouse equation of polymer 
dynamics [|1^], but its behaviour is actually quite different and is in fact extremely unstable 
and would lead to proliferation of length. This can be seen by a transformation to Fourier 
modes = / dsv{s) expjzfcs} gives rfc(t) = rjt(O) exp{A(A;)t — ■iA;(^t} with A(fc) = Ak"^ — Bk^ 
so there is always a band of unstable modes for which A > 0. 

In the form of equation(^) above, the essential dynamics is physically transparent and 
numerical interpretation very efficient. Similar equations (but without noise) such as those 
found in viscous fingering [|2D[ could also be studied in this way. 



We numerically integrated equation(^ above starting from an initial condition of a 
straight line and from Figure 1 it can be seen that meandering patterns are developed. Fur- 
ther evolution of the pattern leads to a plane filling, many time self-intersecting curve |jl3 



We now include the non-local part of the dynamics, the formation of ox-bow lakes; when 
the river intersects itself it loses that part of the curve between the two points of contact. 
This is the way the system is stabilised. We stress that locality here refers to s the position 
along the river and not r the position in space. 



Before we proceed, it is useful to rewrite the EOM as a Fokker-Planck equation [21] so 



that the probability distribution P([r(s)], t) of the points on the river obey the equation 

|p = ap (7) 



where 



O I dL d ^ ( ^ + C + A + B ^^^^ 
Jo 6r(s) \ 6r(s) ds ds^ ds'^ , 



The integral over L acknowledges the fact that the length does not remain constant. Since 
the river is taken as very long, the fact that the length L varies is not so important as we 
assume that near the steady state (L) is a constant in time. 

The oxbow interaction can be considered as a transition between states and we can define 
a transition probability in a Boltzmann type master equation [PTI. The equation will be of 



the form dtPa = OfPa + Z^fe [TabPb — ^6a-Pa] whcrc Ph refers to states scattering to Pa or 
being scattered to from Pa via ox-bow creation and Tab is the 'transition matrix'. For our 
purposes, it is sufficient to make the association P([r(s)],t) = limAr^oo • • • ,i"iv],^) so 
the functional can be taken as the infinite limit of a discrete function of many variables. In 
this notation the creation of an ox-bow happens when the curve intersects itself at i and 
j where without loss of generality, j > i corresponds to the transformation of an initial 
probability distribution Pj for which = rj to a final distribution Pp for which the portion 
of the curve between rj and rj has vanished and so the curve must be re-labelled such that 
k > j, Tfc rfc_(j_j). We can obtain an approximation of the final distribution in terms of 
the initial distribution by a generalisation of a Taylor expansion to functionals so that in the 
continuum notation, the initial and final distributions are given by, respectively, Pj[r{s), L, t] 
and PF[r{s), L — 6L, t + 6t] where 6L is the length of the ox-bow. Therefore the creation of 
an ox-bow by intersection at Si and sj with Si < sj is described by 

Pp[{r(s)}, Si, Sj] - P/[{r(s)}, Si, Sj] ^ 



L . 6Pi / dr 4 d^r ^/ .\\ 



where Sij = \si — Sj\. Now this will only work if we can consider \si — Sj\ to be small and so 
by default this approximation will only be able to deal with small loops, i.e. |sj — Sj| << L. 

The transition probability r will depend on the rate of approach of points r(sj) and r{sj) 
and so r = Xdt{J J dsidsj6[r{si) — r{sj)]} where A is a constant. The sum over b in the 
Fokker-Planck equation corresponds to a configurational sum over the state being scattered 
to or from. 

To determine the effect of ox-bows, we make the assumption that their effect added to our 



prescription will give us, in the statistical sense, a steady state. In a sense, we construct such 
a steady state. The EOM becomes dtPs = OrPs+loxPs where Xox is an effective interaction 
due to the ox-bows. Because we are in the steady state, {dtPs) = so that we may write 
(Ps) oc exp {— F[r]} where F = Fo[r]+F'[r] and the operator acting on P, is given hjIoxPs = 
Jds6/6r{5F'/5vPs) and the other part, 6Fo/6r{s) = +Cdr/ds + Ad^r/ds^ + Bd^r/ds^. 
Fq is the part of the free energy from the local dynamics and we want F' the part from the 
ox-bow to somehow cancel out the instabilities in the local EOM. It turns out that F' has 
two parts. 

The first part of the interaction can be dealt by looking at the average value of the 
transition away from the present state by loss of an ox-bow, and we find an interaction term 

F{ ~ A y dsds'S[r{s) - r(s')] . (9) 

To obtain this we used the fact that {dtr) = SF/Sr{s) and S{Ps)/Sr = 
{Ps)SF/6r and we rewrote the transition probability in the form r = Xdtr{si) ■ 
6/6r{si) {/ / dsidsj6[r{si) — r{sj)]}. This is the same form of term that would be obtained 
from the dynamics of the self-avoiding walk (SAW) Hamiltonian. 

The second part is obtained from the transition into the present state from another by 
loss of an ox-bow. Here we find that the effect of the transition is given by 

5 6F' 5P ( .^.dv .^.d^r \ ,^ . 

where = (A / dsidsjdt{6[r{si) — r(sj)]}|sj — Sj\°') are a simple function of A and mea- 

sures the average rate of approach of points along the curve. The effect of this will be that if 
we consider the original equation we have a renormalisation of parameters due to the inclu- 
sion of the effective field of the ox-bow interaction, A ^ A = A— (72(A) and C ^ C = C~5'i('^)- 
For the system to remain in a stable state, we musthsive A— (72(A) = —2AD/B and ( = (71(A). 
We get the same results from independent scaling arguments. Combining both parts, we ob- 
tain a stable state that should correspond to a system always near that of a SAW confirming 
recent lattice simulations and measurements from field data [E^l. 



We also included the creation of ox-bows to the numeric analysis of the meander dynam- 
ics by a relabelling procedure on contact, and we generated stable meander configurations 
in contrast to the situation when they are not present. In Figure 2 we have the configura- 
tions after 200, 000 time-steps for the system with and without the inclusion of the ox-bow 
mechanism starting from an initial configuration of a straight line in both cases. The mean 



spread of the river path Yg = y {{Y — F)^) as a function of length L was also calculated and 
we obtained the relationship Yg ~ L for short distances and a cross-over to Yg ~ fo^ 



long distance, agreeing with field data ||23 



We now briefly consider some extensions to the model to make it more realistic; these 



will be discussed fully elsewhere . In the expression for the fluid velocity for the channel, 
higher order terms are of the type 6u ~ Kexp{— 7(5 — s')} where 7 is a constant depending 
on the parameters of the flow and infract it can be shown that such corrections lead 



to an addition of a constant vector to the EOM. In other words it gives the river a well 
defined direction and does not change the dynamics described above. It is also instructive 
to consider higher order terms in the Taylor expansion undertaken to describe the ox-bow 
creation. These terms lead to a renormalisation |TB| of the noise so that D ^ D. Finally we 
consider the effect of 'valley confinement' and slope to the river channel. These will result 
in a potential U[r] of the form {g, 0) + (0, hy) added to the EOM. 

In conclusion, we have developed a phenomenological model of river meanders which 
reproduces all the major features seen in nature motivated by an analysis of flow in a 
meandering channel. An essential feature of the model is the competition of ox-bow creation 
with a curvature instability of the river channel. It is particularly interesting that with the 
inclusion of ox-bows to the motion of the river, stable meanders form for a// values of A, B, D. 
The ratio B/A gives the lower length-scale of the river pattern or the size of the smallest 
ox-bows and the value of D gives the 'time scale' over which the system retains memory of its 
previous behaviour. This explains the visual scale invariance of meandering river systems. 
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FIGURES 

FIG. 1. The evolution of a straight line at t = 0, 4 x 10^, 6 x 10^, 8 x 10^ and 10^ (bold line) 
iterations with A = 0.005, B = 0.000012 and D = 0.002 via equation(5). 

FIG. 2. Configuration after 200,000 iterations (a) with oxbows and (b) without oxbow creation. 
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